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Abstract: We consider the mass measurement at hadron colliders for a decay chain of 
two steps, which ends with a missing particle. Such a topology appears as a subprocess 
of signal events of many new physics models which contain a dark matter candidate. 
From the two visible particles coming from the decay chain, only one invariant mass 
combination can be formed and hence it is naively expected that the masses of the three 
invisible particles in the decay chain cannot be determined from a single end point of the 
invariant mass distribution. We show that the event distribution in the \og(E\T / 'Ezr) 
vs. invariant mass-squared plane, where E\t, E%t are the transverse energies of the two 
visible particles, contains the information of all three invisible particle masses and allows 
them to be extracted individually. The experimental smearing and combinatorial issues 
pose challenges to the mass measurements. However, in many cases the three invisible 
particle masses in the decay chain can be determined with reasonable accuracies. 
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1 Introduction 

Two major questions that particle physics are facing today are the origin of the elec- 
troweak symmetry breaking and the identity of the dark matter in the universe. The 
answers to both questions may lie on the physics at the TeV scale. For the first question 
the TeV scale is the scale of the electroweak symmetry breaking. It is widely expected 
that there will be new states at the TeV scale to stabilize the electroweak scale, solv- 
ing the hierarchy problem of the Higgs sector. The second question is related to the 
TeV scale because the most promising dark matter candidate is a new stable weakly 
interacting massive particle (WIMP) at the weak scale. The thermal relic left from the 
Big Bang for such a particle provides the right amount of dark matter in the current 
universe. It is quite possible that these two questions are related and have their com- 
mon origin from the TeV scale physics. The Large Hadron Collider (LHC) is currently 
running and searching for new physics at the TeV scale intensively. Discoveries may be 
made at any time to provide us important clues to these fundamental questions. 

Even though we do not know the exact new theory at the TeV scale, a typical 
collider signature for new physics containing the dark matter is the missing transverse 
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energy {fir) in an event. To ensure the stability of the dark matter particle, a new 
symmetry is often needed such that the dark matter particle is the lightest particle 
charged under this new symmetry while all standard model (SM) particles are neutral. 
There can be other new particles at the TeV scale which are charged under this new 
symmetry, as it is often the case in a more complete theory, including the most popular 
supersymmetric standard model (SSM) where the i?-parity plays the role of the the 
new symmetry [1-4]. There are many other models which share the same feature, e.g., 
universal extra dimensions (UEDs) [5-7] and little Higgs theories with T-parity [8- 
10]. These new particles are necessarily pair-produced at the collider due to the new 
symmetry. After being produced, they decay down to the dark matter particles which 
escape the detector, yielding missing energies, together with SM particles coming from 
the decays. The standard signatures for this class of models are therefore jets and/or 
leptons with 

Missing transverse energy is an important channel to search for new physics at the 
colliders. Events with large missing transverse energies are easy to identify. However, 
once such signals beyond the SM backgrounds are discovered, it is a non-trivial task 
to identify what these new physics events are and to reconstruct their kinematics. 
With at least two missing particles in each event, in general one cannot reconstruct 
the full kinematics on an event-by-event basis without additional information. Some 
more sophisticated techniques need to be developed if one wants to determine the 
properties of the invisible particles including the dark matter particle which appear in 
the process. One of the first things that we would like to know about the new particles 
are their masses. Once the masses are determined, they can be used to reconstruct 
the kinematics of the events, and hence to help determining other properties of these 
new particles. In the past decade, there have been many techniques developed to 
determine the masses of the new particles in the decay chains which end up with 
missing particles. (See Ref. [11] for a review.) The starting point is to isolate the 
signal events that predominantly come from a particular topology 1 and then find the 
appropriate kinematic variables or constraints which depend on the masses of the new 
particles. Some variables, e.g., the invariant mass end point [13, 14] and the stransverse 
Mt2 [15-17]) can be used in a wide variety of event topologies, but in general the best 
mass determination may require different techniques for different event topologies. 

In recent years many studies have been focussed on symmetric decay chains [18- 
31]. The advantage of considering symmetric decay chains is that there are additional 
kinematic constraints coming from mass equalities between the particles in the two 

The event topology can be tested by looking at distributions of various kinematic variables such 
as the invariant masses [12]. 
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Figure 1. Decay chain 



decay chains and the total missing transverse momentum, which allowed accurate mass 
determination even for relatively short decay chains. On the other hand, in a model 
where many different decay patterns can happen, there are often many more events 
with asymmetric decay chains. There have been attempts to generalize the techniques 
developed for symmetric decay chains to asymmetric chains, but so far only to different 
mother or daughter particle masses between the two chains, and some prior knowledge 
about the asymmetry of the events is required [32, 33]. To be complete general, it is 
worthwhile to find kinematic techniques or variables which can be applied to a single 
decay chain. Then they can be used even in a collection of different types of signal 
events as long as they contain one identical decay chain. In the case of limited statistics, 
they may give the first estimate of the masses before there are enough symmetric-chain 
events for analysis. 

In fact, the early attempts of the mass determination in SUSY-like events were 
based on single decay chains [34-40]. However, to obtain enough kinematic constraints, 
long decay chains with at least 3 steps are required. There are many invariant mass 
combinations which can be formed with 3 or more visible particles, which provide 
enough mass relations among the invisible particles in the decay chain. So far there is no 
technique for extracting the masses in a shorter decay chain with only two or less visible 
particles. Naively one may imagine that this task is impossible. For example, consider 
a decay chain starting with a mother particle Y which decays through an intermediate 
on-shell state X and ends up with the missing particle N as shown in Fig. 1. There is 
one visible particle from each step of the decays and the visible particles are labeled as 
1 and 2 in the figure. There are three invisible particles (Y, X, N) in the decay chain 
but there is only one invariant mass that can be formed from the two visible particles. 
One invariant end point certainly cannot determine three masses, but only provides 
one relation among them. It is also well known that the shape of the invariant mass 
distribution depends on the spin of the intermediate particle X and the chiralness of its 
couplings, but not the masses [41-53]. One needs to find two more independent mass 
relations in order to solve for the three masses. In this work we make the first attempt 
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in a model-independent way to measure the invisible particle masses of the short single 
decay chain of Fig. 1. We show that in certain cases, it is possible to determine all 
three masses in the decay chain. We are not aiming for high accuracies. After all, this 
is a difficult case so even a 30%-50% measurement is infinitely better than not being 
able to determine them at all. Any information obtained or techniques developed here 
can also be used in more complicated event topologies which contain this decay chain 
as a subprocess to provide additional constraints. Furthermore, the variables used for 
mass determinations are often useful to separate signals from backgrounds and hence 
can be used in new physics search in the first place. 

A crucial observation is that the masses are mostly determined by certain special 
events (and their nearby events). Those events often lie at the end point or the peak of 
some kinematic distribution. Mathematically, those special events make the kinematic 
constraint equations degenerate. In order for the degenerate equations to have solu- 
tions (since the correct masses should be compatible with those events), this implies 
certain relations among the coefficients of the constraint equations, which depend on 
the invisible particle masses. In this way we obtain mass relations among the invisible 
particles. It has been shown that the common variables such as the end points of the 
invariant mass and the transverse mass variables M?, M?2 can all be understood in 
this way [54]. There some constraint equations used are quadratic so the special events 
lie at the end point of a kinematic distribution. The regions near these points are often 
where signal events accumulate due to the phase space restriction, so sometimes they 
can also be used to select signals over backgrounds. 

For the decay chain shown in Fig. 1, the constraint equations are 

Py = m Y , (1.1) 
(PY-P2) 2 =m 2 x , (1.2) 
(p Y -p 2 -Pi) 2 = m 2 N , (1.3) 

where p 2 and p\ are the 4-momenta of the visible particles from the decays of Y and X 
respectively and for simplicity we take them to be massless p\ = p\ = which is a good 
approximation for most SM visible particles. Since we do not use the information from 
the other decay chain, there is no constraint from the missing transverse momentum. 
Taking the differences of these three equations, we can obtain two linear equations in 
the unknown momentum py , 

2p 2 py = m Y - m\ = A 2 , (1.4) 
2p x p Y - 2p 1 p 2 = m\ - m 2 N = A 1; (1.5) 

where A lj2 are defined as the corresponding mass-squared differences. Together with 
any one of the quadratic equation above, we have all the independent kinematic con- 
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straints. From the discussion in the previous paragraph we should look at the events 
which make these equations degenerate. (For the quadratic equation we just take the 
tangent on the curve.) The events with visible momenta pi, p 2 which make the three 
equations (including the tangent of the quadratic equation) degenerate can be shown 
to satisfy 

{p 1+ p 2 f = «- m x)(fx-m%) = AiA, 

m x m x 

which is nothing but the well-known end point of the invariant mass distribution of the 
visible particle system, m\ 2ms ^ = (pi +j>2) max - By locating the end point of the invari- 
ant mass distribution, we obtain one relation among the three masses my, mx, tun, 
eq. (1.6). 

Because this decay chain can have any boost and orientation in the laboratory 
(lab) frame, one might think that the Lorentz-invariant mass is the only meaningful 
quantity in this process and eq. (1.6) is the only mass- dependent kinematic variable 
that we can get. However, by examining the two linear equations, eqs. (1.4), (1.5), one 
can see that there is another case where the two equations become degenerate, that 
is, when the two visible particles are parallel, p± oc p 2 . In this case the invariant mass 
vanishes, (pi + p 2 ) 2 = 2pip 2 = 0. If we take the ratio of the two linear equations in 
this case, we have 

where E\ and E 2 are the energies of the particle 1 and particle 2 respectively. Because 
the two 4-momenta are parallel, the ratio is invariant under any Lorentz transformation, 
and it gives the second relation among the three unknown masses. 2 One may question 
whether such a co-linear degeneracy is useful in practice, because it rarely happens and 
depending on the nature of the visible particles, the two visible particles may not be 
resolvable in a real experiment. So, instead of using only the extremely rare special 
events when the two visible particles are parallel, we will look at the event distribution 
in the E\jE 2 vs. m\ 2 space (more precisely Eyr/E-xr where the subscript denotes the 
transverse energy as it is invariant under longitudinal boosts). We find that the two- 
dimensional distribution contains useful information and under certain circumstances, 
it can even provide one additional mass relation which allows us to solve for all three 
invisible particle masses. Even in the cases where three masses can not be independently 
and accurately determined, the two-dimensional event distribution provides information 
about the decay chain beyond what is contained in the one- dimensional invariant mass 
distribution. 



2 This relation was also used in Ref. [55]. We thank M. Nojiri for bringing it to our attention. 
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Y 


X 


N 


2 


1 


particle 


left-handed down squark 


2nd chargino 


anti-sneutrino 


up quark 


electron 


mass[GeV] 


777 


465 


292 









Table 1. A summary of the decay chain studied in Section 2. The visible Standard Model 
particles are treated as being massless. 



This paper is organized as follows. In section 2 we analyze at the parton level the 
event distribution of the decay chain in the \og(E\T / E2T) vs. invariant mass-squared 
plane. We show that the masses of the three invisible particles can be obtained by 
fitting the distribution with a simple curve. Some more examples are presented in 
appendix A to show that the method applies to a wide range of models and exceptions 
are included in appendix B. In section 3, we consider the method in more realistic 
situations by including experimental smearing effects and combinatorial backgrounds. 
The case when the order of the two visible particles cannot be identified on an event- 
by-event basis is discussed in section 3.1. The combinatorial problems in the case when 
there are other visible particles identical to one of the visible particles from the decay 
chain in the same event are discussed in section 3.2. Conclusions are drawn in section 4. 

2 The Kinematic Distribution for the Single Decay Chain 

In the Introduction, we have argued that for the decay chain in Fig. 1, in addition to 
the invariant mass end point, a second mass relation can be obtained by looking at the 
energy ratio of the events where the two visible particles are parallel. However, such 
events are rare and may not be usable in real experiments. To avoid these problems, a 
natural thought is to look at the "nearby" events to see if one can extrapolate to the 
point that we are interested in. To do that we should examine the event distribution 
in the two-dimensional space of the energy ratio of the two visible particles and their 
invariant mass. 

To study this with a concrete example, a sample of events from such a decay chain 
are generated with the SUSY LM2 point [56, 57] chosen as the underlying model. The 
masses of the particles in the decay chain are given in Table l. 3 10 4 events are generated 
using MadGraph 4.4.49 [61] at the parton level to reduce the statistical fluctuations. In 
this section we do not include experimental smearing and backgrounds, and assume no 
combinatorial problem. These issues will be discussed in the next section when we deal 

3 This point has been ruled out at the LHC [58-60]. However, we just use this point for the purpose 
of illustration. The method that we develop is independent of the overall mass scale of the spectrum. 
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Figure 2. The histogram of the invariant mass-squared of the two visible particles. In the 
massless limit (pi +P2) 2 = %PiP2- The vertical line indicates the value of ^f' 2 . 
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Figure 3. The scatter plot of log (Ei/E 2 ) vs. 2p±p 2 for all the events. 



with realistic experimental situations. The invariant mass-squared distribution of the 
two visible particles is shown in Fig. 2. The distribution has an end point which can 
be clearly identified. The triangular shape of the distribution is a characteristic of the 
spin-1/2 intermediate state (chargino) in the decay chain. Now we would like to look 
at the distribution of the energy ratios. We take the logarithm of the energy ratio to 
make it more symmetric between the two particles. The two-dimensional distribution 
in the space of log(Ei/E 2 ) vs. invariant mass-squared 2pi p 2 is shown in Fig. 3. One can 
clearly see some interesting pattern of the distribution. As expected, the distribution 
of log (Ei/ 'E 2 ) converges to a point at log(Ai/A2) when the invariant mass goes to 
zero. Away from that point, the \og(Ei/E 2 ) distribution spread out and for a fixed 
invariant mass, it is more or less symmetric about some center point which moves up as 
the invariant mass increases. This qualitative feature already tells us some important 
information about the decay chain. It allows us to figure out which visible particle 
comes from the first step decay and which comes from the second in the case where the 
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two visible particles are distinct, because if we switch E\ and E 2 the distribution will 
move down instead. In our example it means that the quark jet being emitted before 
the lepton can be determined rather than assumed. 

To understand this distribution, let us rewrite eqs. (1.4), (1.5) as 



ZPiPy 

2P2PY 



Ai + 2p 1 p 2 , 
A 2 , 



(2.1) 
(2.2) 



and take the ratio between these two equations. We have 



Ai + 2pip 2 


_ PiPy 


_ E x Ey - 


-Pi 


Py 


_ E x Ey 


- \Pi\\Py\ cos6 1Y 


A 2 


P2PY 


E 2 Ey ' 


-P2 


Py 


E 2 Ey 


- \pi\\Py\ cos 2Y 



lab 



ExEy{1 - Py cos6 1Y ) 



E 2 E Y (l-f3 Y cos6 2Y ) 
Ei(l — (3 Y cos Oiy 



lab 



E 2 (l -(3 Y cos9 2Y ] 



(2.3) 



lab 



where the cos6 , i( 2 )y is the angle between particle 1(2) and particle Y, the subscript 
"lab" indicates that the angles are measured in the lab frame and Py is the magnitude 
of the velocity (boost) of particle Y, defined as /3y = \py\/E y . It is easy to see that 
eq. (2.3) reduces to the simple relation, A x /A 2 = Ei/E 2 , when particle 1 & 2 are 
parallel to each other. Now taking the logarithm of eq. (2.3) we obtain 



E x Ai + 2pip 2 l-/3 y cos6» 2 y 

log — = log + log 

i±i 2 i\ 2 

! Ai + 2 PlP2 
= log 7 1- log 



1 — f3y cos 6iy 
1 + (3y cos Oiy 



1 + Py cos 2Y 



lab 



Y 



(2.4) 



where the subscript u Y n denotes that the angles are measured in the rest frame of 
particle Y. More specifically, cos6 l i( 2 )y| y is the angle between particle 1(2) (measured 
in the rest frame of particle Y) and particle Y (measured in the lab frame). From the 
first line to the second line of eq. (2.4) we simply perform a Lorentz transformation 
and use (assuming particle 1 & 2 are massless) 



cos6>i( 2 )y 



cos6>i (2) y + Py 



lab 



1 + Py cos 0i(2)y 



(2.5) 



Y 



The reason for writing the expression in terms of the angles in the rest frame of particle 
Y is that cos^^ly is directly related to the polarization of particle Y. 

The left-hand side of eq. (2.4) can be measured experimentally. The first term 
on the right-hand side of the equation involves the unknown masses to be determined 
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Figure 4. The distribution of log ^ ^ | ^ as a function of the invariant mass-squared. 

The average of log i+^y cos ely I Y * s ~ O-^ll- It is roughly distributed evenly around zero and 
peaks at zero. 

and the invariant mass of the visible particles which is also measurable. The second 
term on the right-hand side, on the other hand, is not measurable on an event-by-event 
basis. It involves the unknown momentum of the invisible particle Y. When we plot 
the events on the log^i/-^) vs. 2pip 2 plane, this term causes the spread in the vertical 
direction. Nevertheless, if the directions of particles 1 and 2 measured in the rest frame 
of Y are not correlated with the direction of Y itself, we expect that the second term 
will be evenly distributed around zero and peak at zero. Fig. 4 shows the distribution 
of log ! + o y COS a iy L as a function of the invariant mass-squared and one can see that 

° 1+Py cos #21" \ Y ^ 

this is roughly true for any invariant mass. In this case, we can fit the distribution with 
a two-parameter curve 



which has the least total \ 2 measured in vertical distances. Consequently, the two 
fitted parameters A 1; A 2 give estimates of the two mass-squared differences A 1; A 2 . 
As a result, A 1 and A 2 can be determined individually, not just their ratio. Combined 
with the end point of the invariant mass-squared, AiA 2 /m^, they can be inverted to 
solve all 3 unknown masses. 

Actually, eq. (2.4) can be applied in any longitudinal frame as we can give this 
system an arbitrary longitudinal boost. In particular, one can boost each event to a 
frame where E1/E2 = E\t/E2t- This in general removes a large longitudinal boost of 
the mother particle Y which happens when the center of mass of the collision is highly 
boosted. Therefore, we can use the distribution in the \og(En< / E 2 t) vs. 2pip 2 space 



log — = log 



2 



Ai + 2 Pl p 2 
A 2 



(2.6) 
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Figure 5. The scatter plot of log (E\t/E2t) vs. 2pip2 for all the events. The red curve 
is log (E\t/E2t) = log ( Al "^J lP2 )■ The blue curve is a least square fit to the data with the 
same function, treating Ai and A2 as unknown parameters. The standard deviations for all 
the points are assumed to be the same. 



instead. Such a distribution is invariant under the longitudinal boosts and hence it is 
less sensitive to the center of mass of the production mechanism. It is also obvious 
for parallel massless visible particles, E\t/E2t = E1/E2. As shown at the end of this 
section, we find that the fit in the transverse energy ratio space works somewhat better 
than in the total energy ratio space for most of the cases, so we will use the transverse 
energies in our variable log(EiT / E 2 t) throughout the paper. 

Fig. 5 shows the fitted result of the scatter plot in the \og(Ei T / E 2 t) vs. 2p x p 2 plane. 
The red curve is the function 

log tt- = log r (2.7) 

&2T ^2 

for true Ai and A2. The blue curve is the least square fit to the data of the same 
function (2.7) but treating Ai, A2 as fitting parameters. We can see that indeed the 
two curves are very close to each other. The term log IXfiy cos bIy I y * s l nc l eec l evenly 
distributed around zero in this frame for this case. Assuming that we can measure the 
invariant mass end point accurately, we can then reconstruct the 3 invisible particle 
masses my, rax, rax with the fitted values for Ai and A 2 . The result is shown in 
Table 2. We see that, compared to the true values, the reconstructed masses are quite 
accurate. 

We have checked this mass reconstruction method for models with different mass 
spectra and different particle spins in the decay chain, and we find that it works well 
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Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log(Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


true 


1.310 x 10 5 


3.875 x 10 5 


-1.08 


777 


465 


292 


reconstructed 


1.370 x 10 5 


3.838 x 10 5 


-1.03 


780 


473 


295 


error 


+4.6% 


-0.96% 


+5.5% 


+0.34% 


+1.8% 


+1.0% 



Table 2. True values, reconstructed values and the errors of the six quantities for the fit 
to all the events. The errors are calculated using roconstm^d-truo ( except for i og (A 1 /A 2 ), 
which is reconstructed — true) and do not represent the statistical fluctuation. In solving for 
the masses, we have used the true invariant mass end point value assuming that it can be 
accurately determined. The uncertainly in determining will add additional error on the 
reconstructed values. 



for a wide range of different models and spectra. A few more examples can be found 
in Appendix A. However, there are two cases where our method falters: 

1. The mother particle Y is polarized and preferentially emits particle 2 in the 
forward or backward direction when it decays. In this case it is clear that cos# 2 y | y will 
have a nonzero average value. On the other hand, the direction of the particle 1 is less 
correlated with the polarization of the mother particle Y. The term log \~^y cos ely I y 
will no longer be distributed evenly around zero, but has a bias depending on the 
cos# 2 y| y distribution. By Taylor expanding the expression 

Ai + 2 PlP2 Ai 2p lP2 

lQg A 2 = bg A- 2 + "AT + ' ' ' ' (2 ' 8) 

we see that the curve has the intercept log(Ai/A 2 ) at 2 Pl p 2 = and the slope 1/Ai to 
the first order. If (cos# 2y L) > 0, we will obtain a fitted Ai (also A 2 since their ratio 
is fixed by the intercept) greater than the true value, then the reconstructed masses 
will be too large, and vice versa. 

2. The mass difference between two invisible particles (in particular between Y 
and X) is small. The visible particle coming from the decay will be soft. There are 
two effects which can affect the mass determination. First, if the mass difference is not 
much larger than the width of one of the particles, the on-shell approximation does not 
work well. The off-shell effects can be asymmetric, which makes the "effective" value 
of A 2 or Ai significantly different from the true value and introduces a bias to the term 
l°g ( Al + ^ lP2 ) i n e Q- (2-4). Second, if A 2 is small and the typical energy of the visible 
particle 2 is not much larger than the Et cut (which is necessary in a real experiment), 
those emitted in the forward direction of particle Y have a larger chance to pass the the 
trigger and the Ex cut than those emitted in the backward direction. This introduces a 
"fake" forward polarization to particle Y, which also causes a problem as we discussed 
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Figure 6. Histograms of log {E\tI^it) for different ranges of invariant masses. The chart 
on the righthand side indicates the range of invariant masses as a multiple of for each 

histogram. The histograms have approximate Gaussian shapes. 

in point 1. In practice, the mass determination is more of a problem for a small A2 
than a small A^ 

In these cases, our simple method does not give an accurate determination of the 
three invisible particle masses. However, the ratio A1/A2 is still usually well determined 
by the interception of the fitted curve with the axis of the zero invariant mass. If a 
third mass relation can be obtained in some other ways (e.g., double-chain events), this 
distribution can provide some other non-trivial information about this process such as 
the polarization of the particle Y. More detailed studies of these two cases will be given 
in Appendix B. 

For the fit in Fig. 5, every event contributes with the same weight. In reality, 
backgrounds and noises are always present. There is no reason to expect that the 
background events will be distributed evenly around the log ( AlH ^ iP2 ) curve. The 



backgrounds can hence cause a bias in fitting. On the other hand, as long as background 
events are subdominant in the event sample, the central peak of the log(£'iT/-E l 2T) 
distribution for a fixed 2pip 2 may not be significantly affected. Therefore, a curve 
fitted to the peak locations of the log(-Err / ' E 2 t) distributions for different invariant 
masses may be less sensitive to backgrounds than a direct \ 2 n ^ to all events. We 
would like to check whether such a procedure gives the result as good as the result in 
Table 2 in the absence of backgrounds first. 

In Fig. 6 we divide the events into five sets with different ranges of invariant masses. 
For each set of events, the histogram of log (E±t / E 2 t) has an approximate Gaussian 
shape. This indicates that we could divide the events according to their invariant 
masses and for each set we extract the peak value using a Gaussian fit. After obtaining 
the peak point for each set, we then fit these peak points to the curve log (Ei T /E 2T ) = 
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Figure 7. The events are divide into 20 sets with equal width (0.05 x ) of invariant 
mass. For each set we extract the peak value by fitting with a Gaussian distribution. The 
horizontal coordinate (2p\p2) of each point is the middle point of each division. By doing this 
we have 20 points of peak values and we fit it with the curve log (E\t/E2t) = log ( AlH ^ lP2 ), 
treating Ai and A2 as unknown parameters. The peak points and fitted curve are shown in 
blue. The red curve is log {EytJE^t) = log ( Al +^ lP2 ) with true values of Ai and A2. The 
error bar of each points is estimated by the formula a[x] = ^j=, where a[x] is obtained from 
the Gaussian fit and N is the number of events in the set. 





Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log (Ax/As) 


m y [GeV] 


m x [GeV] 


m N [GeV] 


true 


1.310 x TO 5 


3.875 x 10 5 


-1.08 


777 


465 


292 


reconstructed 


1.343 x 10 5 


3.827 x 10 5 


-1.05 


776 


468 


291 


error 


+2.5% 


-1.2% 


+3.7% 


-0.17% 


+0.64% 


-0.33% 



Table 3. True values, reconstructed values and the errors of the six quantities for the fit 
to the peak points. The errors are calculated using reconstr "ctcd-true ( except f or l g(Ai/A 2 ), 
which is reconstructed — true) and do not represent the statistical fluctuation. 



l°g ( Ai "a^ iP2 )• The result of this fit is shown in Fig. 7, where the events are divided into 
20 sets with equal width (= 0.05 x AiA2/m^) in invariant mass-squareds. The error bar 
of each point is estimated by the formula a[x] = a[x]/\/~N, where a[x] is obtained from 
the Gaussian fit and N is the number of events in that set. The reconstructed masses 
obtained from this fit is shown in Table 3. We see that the results are comparable to 
the ones in Table 2. 

The errors in Table 2 and Table 3 could be a combination of both systematic and 
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l ° mlE2) log(£,/£ 2 ) 




Figure 8. The scatter plot (left) and the peak-point plot (right) of log (E1/E2) vs. 2p\p 2 
for the same set of events as in Fig. 5 & Fig. 7. In both case the red curve is log (E1/E2) = 
l°g ( Al+ ^ lP2 ) with true Ai and A2 and the blue curve is the fitted one. 





m y [GeV] 


m x [GeV] 


m N [GeV] 


true 


777 


465 


292 


using transverse energies 


reconstructed 


764 ± 11 


456 ± 12 


280 ±11 


error 


-1.6% ±1.5% 


-2.0% ± 2.6% 


-4.2% ± 3.9% 


using actual energies 


reconstructed 


821 ±33 


511 ±34 


336 ± 32 


error 


±5.7% ± 4.2% 


±10.0% ± 7.4% 


±15.1% ±11.1% 



Table 4. The reconstructed masses from 5 sets of events, each with 10 4 events, in the 
form of mean ± standard deviation. The reconstruction is obtained by fitting with the peak 
points, using both the transverse energies and the actual energies. The standard deviations 
are estimated using the unbiased estimator (with a factor of J jfz^) throughout this pa- 
per. The results suggest that using transverse energies results in smaller fluctuations in the 
reconstructed masses. 

statistical errors. Nevertheless, they are quite small and both reconstructions are good. 
To further verify the goodness of our results, we repeated the event generation 5 times 
(with 10 4 events for each set), and applied fits to the the peak points for the 5 sets of 
events. The mean and the standard deviation values of the reconstructed masses for 
the 5 sets are shown in Table 4. We also compare with the results obtained from the 
same procedure on the same 5 sets of events, but using the ratio of the actual energies, 
log(Ei/ E 2 ) in the lab frame, instead of the transverse energies. Fig. 8 shows the scatter 
plot and the fit to the peak points of log (Ei/E 2 ) vs. 2pip 2 for the same set of events 
as in Fig. 5 & Fig. 7. The results from the 5 sets of events using the actual energies 
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are also listed in Table 4. We can see that using transverse energies gives better results 
and smaller fluctuations in the reconstructed masses. 

3 More Realistic Mass Measurements at Colliders 

In Sec. 2 we saw that the event distribution in the transverse energy ratio vs. invariant 
mass-squared space can be used to determine all the invisible particle masses in the 
decay chain. It works very well under quite general conditions at the parton level. 
However, there are many complications in performing such a measurement in a real 
experiment. First, we must have a relatively clean sample of signal events to start 
with. Selecting signal events from backgrounds is non-trivial and depends on the type 
of signal events. This is beyond the scope of this paper and we assume that it can 
be achieved for the cases that we are interested in. Even if we have a pure sample 
of signal events, in general we have to face some combinatorial problems. If the two 
visible particles are not distinct or can appear in either order in the decay chain, e.g., 
two leptons from a heavier neutralino decaying down to a lighter neutralino through 
a slepton in a SUSY theory, we do not know which order to take in the energy ratio 
\og(EiT / E 2 t) ■ Another combinatorial issue is that there can be other experimentally 
identical particles appearing in the signal events and there is no absolute way to select 
the correct one. For example, if one of the visible particle is a jet, then we do expect 
other jets to be present in the same event, which can come from the other decay chain or 
initial state radiation (ISR). We will discuss these two types of combinatorial problems 
in the following two subsections. In subsection 3.1, we consider an example of signal 
events with 2 leptons of the same flavor and opposite charges, which can appear in 
either order. We will see that mass determination still works for some cases but not 
for the others, depending on the mass parameters. In subsection 3.2, we consider a 
decay chain which produces 1 jet and 1 lepton with a definite order, but there are 
other jets coming from the other decay chain and/or from the initial state radiation in 
the same event. It is possible that both types of combinatorial problems are present 
simultaneously, as in the case when both visible particles in the decay chain are jets. 
Such cases will be difficult and we do not expect to achieve good mass measurements 
there. 

Finally, experimental smearing of the visible particles from detector resolutions, 
fragmentation and hadronization in the case of a jet, will also deteriorate any mass 
measurement. The smearing effect is more important for jets than leptons. To take 
into account the experimental smearing effect, the parton level events are smeared 
according to the Gaussian errors listed in Table 5 in our studies in this section. They 
roughly correspond to the performance of the CMS detector [62-64]. 



- 15 - 



leptons: 




M < 2.5, p T > 10 GeV, 




^ = 0.008 © 0.00015p T , 




59 = 0.001, 5<p = 0.001. 


jets: 




\rj\ < 5.0, p T > 20 GeV, 


8E T _ 

Srj 


©^#©0.033, for M<1.4, 
IK ©^f© 0.043, for M>1.4, 
= 0.03, 50 = 0.02 for \rj\ < 1.4, 


5rj 


= 0.02, 50 = 0.01 for \rj\ > 1.4. 



Table 5. Parton level events are smeared according to the above Gaussian errors. The cuts 
on px and r\ are consistent with the default cuts in MadGraph. The observables of energy 
dimension are in GeV units and the angular and the rapidity variables are in radians. 

3.1 Combinatorial problems of the 2-lepton signal 

In this subsection we consider the combinatorial problem between the two visible par- 
ticles. As jets will have other problems, we consider that the two visible particles in the 
decay chain are leptons of the opposite charges and the same flavor. This happens in 
SUSY if a heavier neutralino decays through a slepton to the lightest neutralino. The 
two leptons emitted can be in either charge order because the neutralinos are Major ana 
particles and the intermediate state can be either a slepton or an anti-slepton. To focus 
on this combinatorial problem, we assume that there is no other lepton of the same 
flavor in the event. Because we do not know which lepton is particle 1 and which is 
particle 2 on an event-by-event basis, we can not measure \ogiEyr j E 2 t) but only its 
absolute value | \og(E\T I E%t)\ (i-e., the ratio between the larger and the smaller 
E T ). In other words, the scatter plot in the log (E 1t /E 2 t) vs. 2pip 2 plane is folded 
along the log (Eyt/E2t) — axis. This would certainly make the reconstruction more 
challenging. 

Because the scatter plot is folded, the position of the distribution become crucial for 
reconstruction. If the center of distribution is far away from the log (Eit/Ezt) = line, 
the folding will not cause too much of a problem. One can still easily identify the peak 
of the I log(EiT / E%r)\ distribution of each invariant-mass interval, and the slope of the 
fitted curve connecting the peak points tells us the sign of log(Ai/A 2 ). On the other 
hand, if the center of distribution is close to the log (Eit/E^t) = line, the folding 
makes it difficult to identify the pattern and the peak of the unfolded distribution, then 
reconstructing A x and A 2 individually by fitting a curve through the peak points may 
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m y [GeV] 


m x [GeV] 


m7v[GeV] 




case 1 


468 


187 


140.5 


-2.5 


case 2 


468 


237 


140.5 


-1.5 


case 3 (LM2) 


468 


304 


140.5 


-0.56 


case 4 


468 


373 


140.5 


0.40 


case 5 


468 


407 


140.5 


1.0 



Table 6. A summary of the mass spectra in subsection 3.1. The leptons (particles 1 & 2) 
are treated as massless. 

become impossible. 

As the event distribution is folded, it is not good to fit a curve with minimal \ 2 on 
the distribution directly. We will use the second technique discussed in the previous 
section by dividing the data into small intervals of invariant mass-squareds and then 
finding the peak point of the distribution for each set. To extract the peak point of the 
unfolded distribution, we fit each set with a folded Gassian disribution in log (E\t / ' E 2 t) , 
i.e., the Gaussian distribution also folded along the log {Eyt / 'E 2 t) = point. However, 
it is important to notice that the error for a Gaussian fit, a[x] = ^j=, is no longer 
a good estimate of the uncertainty if the distribution is folded. If the center of the 
distribution is close to \og(E 1T /E 2T ) = 0, the folded distribution may not be sensitive 
to the original peak position at all, and the uncertainty is much larger than that of 
the corresponding unfolded Gaussian distribution. To estimate this uncertainty, we 
first use the maximum likelihood method to fit the distribution to a folded Gaussian 
distribution with two parameters, the mean \i and the standard deviation o of the 
unfolded Gaussian distribution. We then find the contour logL = logL max — 1/2 
in the \x-o plane, and the tangent of this contour parallel to the a-axis corresponds 
(approximately) to the boundary of the 68.3% central confidence interval of \x. This 
method has some limitations, however, as the likelihood method requires knowledge of 
the full model up to a few free parameters. We do not know the exact shape of the 
original distribution, but simply treat it as approximately Gaussian. The estimations of 
the peak point and its uncertainty may be off if the original distribution is not close to a 
Gaussian. Once we obtain the estimations of the peak point and the uncertainty for each 
bin, we can fit the peak points with the folded curve | log [Err j E 2 t)\ = | log ( Al ^ lP2 )| 
to extract the values of Ai and A 2 . 4 

To examine how the mass determination depends on the value of logA 1 /A 2 , we 

4 To account for the asymmetric uncertainties but without too many complications, we assume a 
asymmetric Gaussian distribution for the peak point values and use a modified least square fit, i.e., 
the uncertainty used for the fit depends on whether the curve is above or below the point. 
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consider a SUSY process with a heavy neutralino (particle Y) decaying through an 
on-shell slepton (particle X) to the lightest neutralino (particle N), emitting a pair of 
leptons in the decays. The five examples of spectra that we choose to study are listed in 
Table 6. We fix My, Mn and vary Mx to obtain different values for log A1/A2. The case 
3 corresponds to the LM2 point with Y being the 4th neutralino, but the actual model 
is not important. The events are generated through neutralino pair production and 
the polarization of the neutralino is close to zero. Since the neutralino is a Majorana 
particle, by the CP symmetry its decay is symmetric if all final states are included 
anyway. 

Before studying the combinatorial problem we first look at the effect of smearing. 
For case 3, Figs. 9, 10 and 11 show the invariant mass distributions from the two 
visible particles, the scatter plots and the peak points plots in the \og(E 1T /E 2T ) vs. 
2pip 2 plane before and after smearing. We see that the invariant mass end point is less 
sharp after smearing but its existence is still eminent. The value of the end point may 
be obtained with a template fit and we assume that it can be accurately determined. 
Consequently we will use the true end point value in the mass reconstruction as we 
expect much larger uncertainties coming from other quantities. As the smearing effects 
are small for leptons, the distribution in the log (Eit/E 2 t) vs. 2pip 2 plane does not 
change much. The fitted curve also agrees with the true curve quite well after smearing, 
if the ordering of the 2 visible particles were known. 

We now study the folded distributions for the five cases listed in Table 6. We first 
look at the distribution of 10 events for each case to reduce the statistical fluctuations. 
The experimental smearings are included but they do not have a significant effect as 
we see in the above discussion. 

For case 1 (logAi/A2 = —2.5) the scatter plot and the folded peak-point plot are 
shown in Fig. 12. In the scatter plot one can see that a small portion of the events is 
above the log (Eit/E 2 t) = axis, which will be folded. Nevertheless, as the center of 
distribution is far away from the folding line, the folded Gaussian fit works well for all 
bins. The error bars are larger and some touch the folding axis for the last few points 
because the peak points are closer to the folding axis. In this case a fit to all the peak 
points is quite close to the true curve and we can get a good determination of Ai and 
A 2 . 

Fig. 13 shows the plots for case 2 (logA 1 /A 2 = —1.5). Compared to case 1 the 
center of the distribution moves closer to the folding axis. As a result, the error bars 
are larger in the folded Gaussian fits especially for large invariant masses, as shown in 
the right panel of Fig. 13. It is interesting to notice that there are a few points at large 
invariant mass that are far away from the true curve and have relatively small error bars, 
while the true peak points should be close to zero. We check their distributions and it 
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Figure 9. The histogram of the invariant mass-squared of the two visible particles before 
(Left) and after smearing (Right) for case 3. The vertical line indicates the value of ■ 
Both visible particles are leptons. 




Figure 10. The scatter plots in the log {E\x/E2t) vs. 2p\p2 plane before (Left) and after 
smearing (Right) for case 3. In both case the red curve is log (Eyt/E2t) = log ( Al+ ^ lP2 ) 
with true Ai and A2 and the blue curve is the fitted one. 

logCE, r /& r) log(£ lr /-E2 7-) 




Figure 11. The peak-point plots in the log (-Eit/-E2t) vs. 2p\p2 plane before (Left) and 
after smearing (Right) for case 3 (no folding). In both case the red curve is log {E\t/E2t) = 
l°g ( Ai "a^ 1P2 ) with true Ai and A2 and the blue curve is the fitted one. 

turns out that their unfolded distributions are somewhat asymmetric and hence are not 
Gaussian-like. Therefore, the likelihood method assuming a Gaussian distribution did 
not work very well for these points. The fit will deteriorate if these points are included. 



- 19 - 




Figure 12. The scatter plot (Left) and the folded peak-point plot (Right) in the 
log (Eit/E2t) vs. 2pip2 plane for case 1 (log = —2.5). In both case the red curve is 
log (E±t/E2t) = log { Al+ ^ lP2 ) with true Ai and A2 and the blue curve is the fitted one. In 
the folded plot the purple curve is the folded true curve. 
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Figure 13. The scatter plot (Left) and the folded peak-point plot (Right) in the 
log (Eit/E2t) vs. 2pip2 plane for case 2 (log ^ = —1.5). In both case the red curve is 
log (E±t/E2t) = log ( Ai "a^ iP2 ) with true Ai and A2 and the blue curve is the fitted one. In 
the folded plot the purple curve is the folded true curve. 



However, as the nearby points have large errors and touch the folding axis, the center 
of the distribution for this range of invariant mass should be quite close to the folding 
axis. As the points with large error bars are not very useful in the fitting, we choose 
to fit the curve only in the "good region," which contains the first ten points from the 
left, and a reasonable fit can be obtained. 

Fig. 14 shows the plots for case 3 (log Ai/A 2 = —0.56). The center of the distri- 
bution is quite close to the folding line for most range of the invariant mass and the 
fitted points behaves quite badly, as one can see in the plot that most points have very 
large error bars. A naive fit to all points even gives the wrong sign of log A1/A2, which 
may be figured out by a more careful examination of the folded distribution. Anyway, 
this is a bad case for our mass determination method. 



- 20 - 




Figure 14. The scatter plot (Left) and the folded peak-point plot (Right) in the 
log (Eit/E2t) vs. 2pip2 plane for case 3 (log ^ = —0.56). In both case the red curve 
is log (Eit/E2t) = log { Al+ ^ lP2 ) with true Ai and A2 and the blue curve is the fitted one. 
In the folded plot the purple curve is the folded true curve. 




Figure 15. The scatter plot (Left) and the folded peak-point plot (Right) in the 
log (Eit/E2t) vs. 2pip2 plane for case 4 (log ^ = 0.40). In both case the red curve is 
log (E±t/E2t) = log ( Ai "a^ iP2 ) with true Ai and A2 and the blue curve is the fitted one. In 
the folded plot the purple curve is the folded true curve. 



Fig. 15 shows the plots for case 4 (logA 1 /A 2 = 0.40). Now the center of the 
distribution is above the folding axis. The peak points are closer to the folding axis 
for small invariant masses but the corresponding distributions are also less spread. 
Although many points have large error bars as shown in the figure, a curve fit to all 
points seems to give quite good results. 

Fig. 16 shows the plots for case 5 (logA 1 /A 2 = 1.0). The center of distribution 
is quite far away from the folding axis and the effects of folding are unimportant. 
However, in this case A 2 is small and hence particle 2 is quite soft. As discussed in 
Appendix B.2 our reconstruction method starts to produce some bias even without the 
folding. 

The results of mass determination for the five cases are shown in Table 7. For cases 
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Figure 16. The scatter plot (Left) and the folded peak-point plot (Right) in the 
log (Eit/E2t) vs. 2pip2 plane for case 5 (log ^ = 1.0). In both case the red curve is 
log (E±t/E2t) = log ( Ai "a^ iP2 ) with true Ai and A2 and the blue curve is the fitted one. In 
the folded plot the purple curve is the folded true curve. (Note: The origin is not at in the 
peak plot.) 

1, 2, 4 & 5 we also wish to know how well the method works with a smaller sample 
size. It turns out that in terms of statistical fluctuations, for cases 1 and 2 the method 
works reasonably well with 10 3 events. The results are presented in Table 8. 

For cases 4 & 5, some care needs to be taken. As seen in eq. (2.8), the slope 
corresponds to the value of 1/Ai at the first order. When Ai ^> A 2 , the slope is 
small and a small fluctuation in the slope could cause a large fluctuation in the value 
of Ai (also A 2 since their ratio can be well determined) and the mapping is nonlinear. 
As a result, a larger statistics for cases 4 & 5 is needed to get a good measurement. 
Furthermore, a Gaussian-like uncertainty of the reconstructed slope will not result in 
a Gaussian-like uncertainty in the reconstructed Ai when the slope is close to zero 
due to the nonlinear mapping. In particular, the value of Ai goes to infinity when 
the slope goes to zero. This behavior shows up when we look at case 4 with 50 sets 
of 10 3 events. Among the 50 sets of reconstructed masses, 44 of them are relatively 
close to each other but 6 sets has very large reconstructed masses (my > 1 TeV). Case 
5 has a similar behavior, and we found that a larger data sample (2 x 10 3 events) 
is needed for a reasonable reconstruction. Among 25 sets of 2 x 10 3 events, 5 sets 
has my > 1 TeV and the other 20 sets are relatively close to each other. As the 
uncertainty is non-Gaussian, it would be better to quote the central 68.3% confidence 
region instead of mean ± standard deviation. However, it is hard to obtain the form 
of the probability distribution function without generating a large number of sets of 
events. Here we simply present the means and standard deviations of the reconstructed 
masses excluding the bad sets in Table 8. 
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Ax[GcV 2 ] 


A 2 [GeV 2 ] 


log(Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


case 1 


true 


1.523 x 10 4 


1.841 x 10 5 


-2.49 


468 


187 


140.5 


reconstructed 


1.907 x 10 4 


1.974 x 10 5 


-2.34 


494 


217 


167 


error 


+25% 


+7.3% 


+15% 


+5.6% 


+ 16% 


+19% 


case 2 


true 


3.643 x 10 4 


1.629 x 10 5 


-1.50 


468 


237 


140.5 


reconstructed 


5.456 x 10 4 


2.233 x 10 5 


-1.41 


582 


340 


247 


error 


+50% 


+37% 


+8.8% 


+24% 


+43% 


+75% 


case 3 


true 


7.268 x 10 4 


1.266 x 10 5 


-0.555 


468 


304 


140.5 


reconstructed 


2.336 x 10 5 


1.609 x 10 5 


0.37 


734 


614 


379 


error 


+221% 


+27% 


+93% 


+57% 


+ 102% 


+ 170% 


case 4 


true 


1.194 x 10 5 


7.990 x 10 4 


0.402 


468 


373 


140.5 


reconstructed 


9.732 x 10 4 


6.856 x 10 4 


0.35 


407 


312 


0.0 


error 


-18% 


-14% 


-5.1% 


-13% 


-16% 


-100% 


case 5 


true 


1.459 x 10 5 


5.338 x 10 4 


1.01 


468 


407 


140.5 


reconstructed 


1.625 x 10 5 


5.853 x 10 4 


1.02 


511 


450 


200 


error 


+11% 


+9.7% 


+1.6% 


+9.1% 


+11% 


+42% 



Table 7. The results of mass reconstruction from fitting the folded peak points for the 5 cases 
in Section 3.1 with 10 4 events. The errors are calculated using rcconstmctcd~tmc ( except for 
log (A1/A2), which is reconstructed — true) and do not represent the statistical fluctuation. 



3.2 Combinatorial problems of the 1 jet + 1 lepton signal 

In this subsection we study another type of the combinatorial problem when there are 
other particles which can not be distinguished experimentally from one of the particles 
from the decay chain in the event. A typical example is that the two visible particles 
from the decay chain are one jet and one lepton. In general we expect there will be 
other jets present in the same event. They can come from the other decay chain(s) 
in the event or simply the initial and final state radiations. In this case we face the 
problem of not knowing which jet is the correct one to be paired with the lepton in 
the same decay chain. It is possible that there are also other leptons in the same 
event to cause further confusions. To understand the basic issues of this type of the 
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m r [GeV] 


m x [GeV] 


m N [GeV] 


case 1 (10 3 events) 


true 


468 


187 


140.5 


reconstructed 


527±38 


245+34 


195 + 33 


error 


+13% ± 8% 


+31% ± 18% 


+39% ± 23% 


case 2 (10 3 events) 


true 


468 


237 


140.5 


reconstructed 


616 ± 169 


373 ± 160 


280 + 161 


error 


+32% ± 36% 


+57% ± 67% 


+99% + 115% 


case 4 (10 3 events, 44 out of 50 sets have my < 1 TeV) 


true 


468 


373 


140.5 


reconstructed 


438 ± 62 


347 ± 62 


26 + 90 


error 


-6.4% ± 13% 


-6.9% ± 17% 


-81% ± 64% 


case 5 (2 x 10 3 events, 20 out of 25 sets have my < 1 TeV) 


true 


468 


407 


140.5 


reconstructed 


498 + 119 


437+ 119 


151 ± 160 


error 


+6.5% ± 25% 


+7.5% ± 29% 


+7.8% ± 114% 



Table 8. The results of reconstruction of case 1, 2, 4 & 5. The errors are in the form 
"bias+uncertainty" . The means and uncertainties are estimated from a total number of 
5 x 10 4 events, (i.e. For case 1 & 2, we divided the sample into 50 sets, each with 10 3 events, 
reconstructed the masses for each set and obtained the means and standard deviations of the 
reconstructed masses. For case 4 we did the same but we only used the 44 "good" sets out 
of 50 sets. For case 5 we divided the sample into 25 sets, each with 2 x 10 3 events, and only 
used the 20 "good" sets among the 25 sets.) For case 2, only the first ten points (counting 
from the left, see Fig. 13) are used for fitting. 

combinatorial problem we will keep it simple by restricting ourselves to the case where 
only additional jets are present. The analysis should be readily generalizable to more 
complicated situations. 

The jets coming from decays of some heavy particles and the jets coming from 
initial state radiation have very different characteristics. To consider both possibilities, 
we generated the following process using MadGraph. A left-handed down squark and 
its anti-particle are pair produced with the left-handed down squark decaying through 
the same chain studied in Sec. 2 (with the spectrum shown in Table 1), while the anti 
left-handed down squark decays to an anti-down quark and a lightest neutralino. In 
addition, an extra gluon is produced, either through ISR or FSR or elsewhere. The 
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Figure 17. The histogram of the invariant mass-squared of the three possible combinations 

before smearing (left) and after smearing (right). The vertical line indicates the value of 
AiA 2 
"pr- 
events are smeared according to the Gaussian errors listed in Table 5. The signal of 
this process is 3 jets + 1 lepton + Et and we do not know which jet is particle 2. This 
case is, of course, an over simplified version compared with the real processes at the 
LHC, but it should help us understand how to apply our method to the realistic cases. 

Fig. 17 shows the invariant mass-squared distributions for different combinations 
before and after smearing. For convenience we label the anti-down quark from the other 
chain particle 3 and the extra gluon particle 4. The invariant mass-squared distribution 
of the three different combinations are shown in the same plot. In reality we will not be 
able to distinguish them, but only obtain a distribution in which all three are stacked 
together. 5 Fig. 18 shows the scatter plot in the log (Ext/E^) vs. 2pipj plane where 
j = 2, 3, 4 denotes the jet is particle 2, 3 (from the other decay chain) or 4 (extra gluon). 
For better visualization only 1000 events are plotted. For j = 4, the events are mostly 
in the upper left region due to the fact that most gluons are soft. For j = 3 the events 
are more spread on the 2p±pj axis and have less correlation between log (Eit/Ejt) and 
2p l p j . 

Before smearing, the invariant mass-squared distributions have two distinct fea- 
tures. First, for the correct combinations (p±, P2), the distribution has a sharp edge at 
Al t 2 , while the wrong combinations do not have such an edge. Second, for invariant 

m x 

mass below the edge, the correct combination is more likely to have large invariant 
mass than that of the wrong combinations. The shape of the invariant mass-squared 
is model-dependent. In particular, it depends on the spin of the particle X and the 
chiralness of its couplings. If particle X is a scalar, the invariant mass-squared will have 

5 Thc ISR jet may be identified with some probabilities through some sophisticated methods [65, 66], 
but will not be attempted here. 
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Figure 18. The scatter plots in the log (Eit/Ejt) vs. 2p\p2 plane before smearing (left) 
and after smearing (right) where j = 2,3,4 denotes the jet is particle 2, 3 (from the other 
decay chain) or 4 (extra gluon). 1000 events are plotted. No cut (apart from the default cut 
in MadGraph) is applied. 

a flat distribution below the edge, but compared with the distributions of the wrong 
combinations, the two features still hold. On the other hand, if the left-handed down 
squark in the above process is replaced by an up anti-squark, then the quark and the 
lepton in the decay chains will have opposite helicities (in contrast to the same helicity 
in our example). Consequently they tend to go in the same direction and have a small 
invariant mass. In this case there will be no sharp edge at the end point. However, 
given that we need a good measurement of Al ^ 2 from the end point, we will only con- 
sider cases where a sharp edge is present, then generally the second feature also holds. 
Based on these two features we will use a simple and naive way to select the jet-lepton 
pair by choosing the combination with the largest invariant mass below the edge. 

Because jets suffer more experimental smearing than leptons, the edge is more 
washed out compared to the pure leptonic case in the previous subsection. Neverthe- 
less, the existence of an edge in the invariant mass-squared distribution should still be 
identifiable from the right panel of Fig. 17. We again assume that the location of A ^ 2 
can be obtained from a template fit. Any uncertainty for this quantity will add to the 
total uncertainties of the reconstructed masses but we do not expect it to be the main 
source of the final uncertainties. We will not include it in our following discussion. 

To select the correct jet, we first make the following cuts: 

• A pt cut on the jets is imposed. Such a pt cut is inevitable experimentally to 
reduce the QCD backgrounds. It mostly help to remove the ISR jets. However, if 
the pt cut is too high, it may create a fake polarization for the mother particle Y 
as it will favor the events with Y decaying to particle 2 in the forward direction, 
and hence causes a bias in the mass determination. In our example we require 
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jets to have pr > 50 GeV. 



• We require the jet when combined with the lepton to have an invariant mass- 
squared smaller than the position of the edge, A ^Jt 2 ■ This dominantly reduces 
the wrong combinations coming from the jet of the other decay chain. Because 
of smearing a small fraction of the correct combinations will be removed too. 

If no jet survives the above cuts, the event is dropped. If there are more than one jets 
passing both cuts, we select the jet which forms the largest invariant mass together 
with the lepton. 

After selecting the jet in each event, we can proceed as usual by dividing the data 
points according to their invariant mass-squareds and finding the peak location in the 
log(EiT /E 2 t) distribution for each invariant mass-squared bin. The peak-point plot 
after smearing is shown in Fig. 19. The blue points are for all correct combinations and 
the green points are for selected combinations using the method we described above. 
From the blue points we can see the effects of smearing: most points at lower invariant 
massed are shifted slightly upwards while the points with invariant mass close to 
are shifted downwards. The wrong combinations mostly come from the extra gluon at 
low invariant masses and the jet from the other decay chain at high invariant masses. 
We see that the green points at small invariant masses are shifted up significantly due 
to the softness of the extra gluon. At large invariant masses close to the edge, the green 
points are further shifted down compared with the blue points because in this example 
the jet from the other chain on average has a larger energy than the correct jet does in 
the large invariant mass region. 

The wrong combinations make the mass determination more challenging. We see 
that the peak points at both ends of the invariant mass range deviate from the true 
values significantly. Because the distribution of the ISR jets mostly lies at the upper 
left corner and the distribution of the jets from the other decay chain is flat in the 
\og(Eij< I E 2 t) vs. 2pip 2 plane, the effects of the wrong combinations from both sources 
will make the extracted slope smaller than the true value, which will result in too big 
reconstructed masses. To reduce the bias caused by wrong combinations, one may want 
to use only points in the middle part of the invariant mass-squared interval to perform 
the fit. However, the good range for the fit seems to be model-dependent when we check 
models with different spins and spectra. A fixed range in the invariant mass-squared 
interval (e.g., middle 1/3 of points) which works for some models does not work for the 
others. The best strategy that we come up with so far is to have a fixed size of the 
invariant mass-squared range, but keep the upper and lower ends floating to maximize 
the fitted slope. This is motivated by the fact that wrong combinations tend to reduce 
the slope. The range used for fitting should not be too small to avoid a too large slope 
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Figure 19. The peak-point plot of log (E\t/E2t) vs. 2p\p2 after smearing. The events are 

divide into 20 sets with equal width (0.05 x Al ^ 2 ) of invariant mass and there are 20 points 

m x 

in the plot. The red curve is log {E\t/E2t) = log ( 1+ ^ lP2 ) with true Ai and A 2 . The 
blue points are for all correct combinations and the blue curve is the fitted curve with last 
5 points dropped. The green points and curve are for the selected combinations with the 
method described in Section 3.2. For the green curve, point 7 to 13 are used in the fit. 

coming from the statistical fluctuations, especially when the sample size is small. In 
Fig. 19 we fix the length of the fitting range to be 7 points and find that point 7 to 
13 (counting from the left) gives the largest slope after scanning through all possible 
choices. The result obtained from fitting these points is shown in Table 9. 

The fitting method is somewhat heuristic so one would like to know whether the 
errors obtained in Table 9 is typical and how they depend on the size of the event 
sample. To study that we generate 10 5 events which are smeared according to the 
Gaussian errors listed in Table 5. We first divide the events into 10 sets with 10 4 
events in each set, then reconstruct the masses for each set using the above method. 
The results are shown in Table 10. In practice, 10 4 is a very large number in terms of 
number of events which can obtained in real experiments, so we also check how well 
the method works with less events. To do this, we divide the events into 50 sets, each 
containing 2000 events. The results of the mass reconstruction in this case are shown in 
Table 11 using the same method, except that the length of the fitting range is increased 
to 8 points to compensate for the effect of the larger fluctuations due to the smaller 
sample size. 
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Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log(Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


true 


1.310 x 10 5 


3.875 x 10 5 


-1.08 


777 


465 


292 


true combinations 


reconstructed 


1.411 x 10 5 


3.719 x 10 5 


-0.97 


772 


473 


287 


error 


+7.8% 


-4.0% 


+12% 


-0.68% 


+1.7% 


-1.6% 


selected combinations 


reconstructed 


1.727 x 10 5 


4.418 x 10 5 


-0.939 


876 


570 


390 


error 


+32% 


+14% 


+15% 


+13% 


+23% 


+34% 



Table 9. Results of reconstruction from the fit to the peak points after smearing with 
true and selected combinations. The errors are calculated using reconstr "cted-true (except for 
log (A1/A2), which is reconstructed — true) and do not represent the statistical fluctuation. 





m y [GeV] 


m x [GeV] 


mjy[GeV] 


true 


777 


465 


292 


reconstructed 


931 ± 100 


624 ± 103 


445 ± 99 


error 


+20% ± 13% 


+34% ± 22% 


+52% ± 34% 



Table 10. The reconstructed masses from 10 sets of events, each with 10 events, in the form 
of mean ± standard deviation. The reconstruction is obtained by fitting with 7 consecutive 
peak points which maximize the fitted slope. 





m y [GeV] 


m x [GeV] 


m N [GeV] 


true 


777 


465 


292 


reconstructed 


754 ± 157 


442 ± 162 


265 ± 162 


error 


-2.9% ± 20% 


-4.9% ± 35% 


-9.4% ± 55% 



Table 11. The reconstructed masses from 50 sets of events, each with 2000 events, in the form 
of mean ± standard deviation. The reconstruction is obtained by fitting with 8 consecutive 
peak points which maximize the fitted slope. 

The uncertainties of the results in Table 11 are quite substantial, though they 
already represent significant progress as there was no existing method which can deter- 
mine the invisible particle masses for this topology. We have used a rather simple and 
somewhat heuristic method in treating the combinatorial problem and the results may 
be improved if more sophisticated techniques are employed. Some possible improve- 
ments include: 

1. We use a naive method to select the jet to pair with the lepton. We expect that 
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the correct rate can be improved with a better algorithm. The ISR jet may be 
better identified with more sophisticated methods [65, 66]. One can also employ 
a full likelihood method to distinguish the correct jet from the wrong ones using 
more variables. With better selections we may extend the fitting range to reduce 
statistical uncertainties and model dependence. 

2. Our heuristic method in choosing the fitting range is concerted in such a way 
that it can be applied to a wide range of models. If more detailed information 
is available about the process and the impostor particles from other part of the 
event which contains the decay chain, one can design a fitting procedure best 
suit for the specific set of signal events using all available information and hence 
obtain better results. 

4 Conclusions 

In searching for new physics at colliders there are different approaches. Given a spe- 
cific model, one can design cuts and strategies to obtain the maximal reach for that 
particular model, and perform the most likelihood fit to extract model parameters. 
However, they may not be suitable for a different model. As there is a vast number of 
possibilities for new physics at the TeV scale, one can only perform specific searches 
for a limited number of models and there is no reason to expect that they are favored 
over the other models. On the other hand, model-independent global searches compare 
data and standard model predictions in a global sample of high-p^ collision events and 
look for excess in the high-p T tails as signals for new physics [67]. They do not require 
specification of any particular model. The drawback is that new physics signals may 
be subtle and hide in the correlations of observables, then the global searches are not 
effective to uncover them. 

An intermediate approach is to focus on the possible event topologies of new physics 
without completely specifying the underlying model. An example in this direction is the 
simplified model approach for new physics searches [68]. For any given event topology, 
one can look for kinematic variables and relations which exhibit particular features for 
that topology, so that one can use them to distinguish new physics from the standard 
model backgrounds and also extract relevant parameters of new physics. Because many 
different models can produce signal events of the same topologies, the study of any 
particular event topology can apply to a wide range of possible new physics models. 

In this paper we consider the event topology of a 2-step decay chain ending with a 
missing particle. This topology appears in many models which contain a stable or long- 
lived neutral particle that escapes detection. Out of the two visible particles coming 
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from the decay chain, the obvious kinematic variable to look at is the invariant mass 
combination of the 4-momenta of the two visible particles, and one might thought that 
it is the only relevant kinematic variable for this topology. We show that the logarithm 
of the transverse energy ratio of the two visible particles is also a useful variable. The 
event distribution in the \og{Eyr / E-^r) vs. invariant mass-squared space carry useful 
information beyond what is contained in the one- dimensional invariant mass distribu- 
tion. In many cases it allows the extraction of masses of all three invisible particles in 
the decay chain. The shape of the distribution may also be used to distinguish signals 
from backgrounds. 

The method of searching for new useful kinematic variables may be generalized to 
other event topologies. Digging new physics signals out of the tremendous backgrounds 
at hadron colliders is always a challenging task. Finding the most effective kinematic 
variables for new physics signals will improve the search reaches and help us to figure 
out the properties of the underlying theory once it is discovered. 
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A More Parton-level Examples 

In this Appendix we present a few more parton-level examples for our mass determina- 
tion method. It works well for a variety of models with different mass spectra and spins. 
In all examples in this Appendix, the mother particle Y decays isotropically and the 
order of the visible particles are assumed to be known. Anisotropic decays of particle 
Y will cause a bias in mass determination which will be discussed in Appendix B. The 
combinatorial problems of indistinguishable particles or orders are discussed in sec. 3. 

Example 1 

In this example we consider a SUSY process with a heavy neutralino (particle Y) 
decaying through an on-shell slepton (particle X) to the lightest neutralino (particle 
N). The mass spectrum and spin configuration are summarized in Table 12. The 
invariant mass distribution is shown in Fig. 20. The scatter plot and the peak-point 
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plot are shown in Fig. 21. The result of mass reconstruction from fitting the peak 
points is shown in Table. 13. 





Y 


X 


N 


spin 


1/2 





1/2 


mass[GeV] 


468 


304 


140.5 



Table 12. The summary of the mass spectrums and spin configurations of Example 1 in 
Appendix A. 
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Figure 20. The histogram of the invariant mass-squared of the two visible particles for 



Example 1 in Appendix A. In the massless limit (pi + P2) 

indicates the value of Al ^ 2 . 

m x 



2pip2- The vertical line 




Figure 21. The scatter plot (left) and the peak-point plot (right) of log (Eit / '-Est) v s- 2>PiP2 
for Example 1 in Appendix A. In both case the red curve is log (E\t/E2t) = log ( Al+ ^ lP2 ) 
with true Ai and A2 and the blue curve is the fitted one. 
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Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log(Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


true 


7.268 x 10 4 


1.266 x 10 5 


-0.555 


468 


304 


140.5 


reconstructed 


7.517 x 10 4 


1.295 x 10 5 


-0.544 


477 


313 


150.4 


error 


+3.4% 


+2.3% 


+1.1% 


+1.9% 


+2.9% 


+7.1% 



Table 13. The result of mass reconstruction from fitting the peak points for Example 1 
in Appendix A. The errors are calculated using reconstructcd-tmc ( excep t for log(Ai/A 2 ), 
which is reconstructed — true) and do not represent the statistical fluctuation. The statistical 
fluctuation are expected to be at the same order as the example in Section 2. 



Example 2 

This is a simplified model. The interactions are vector-like so that there is no preferred 
direction for the particle Y decay. The mass spectrum and spin configuration of the 
particles in the decay chain are summarized in Table 14. The invariant mass distribution 
is shown in Fig. 22. The scatter plot and the peak-point plot are shown in Fig. 23. The 
result of mass reconstruction from fitting the peak points is shown in Table. 15. 





Y 


X 


N 


spin 


1 


1/2 


1 


mass[GeV] 


1000 


800 


500 



Table 14. The summary of the mass spectrums and spin configurations of Example 2 in 
Appendix A. 



2p\P2 



Figure 22. The histogram of the invariant mass-squared of the two visible particles for 
Example 2 in Appendix A. In the massless limit (p± + P2) 2 = %PiP2- The vertical line 
indicates the value of Al ^ 2 . 
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Figure 23. The scatter plot (left) and the peak-point plot (right) of log (-Eit/^t) vs. 2p±p2 
for Example 2 in Appendix A. In both case the red curve is log (-Eit/^t) = log ( Al +^P lP2 ) 
with true Ai and A2 and the blue curve is the fitted one. 





Ax[GeV 2 ] 


A 2 [GcV 2 ] 


log (Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


true 


3.900 x 10 5 


3.600 x 10 5 


0.080 


1000 


800 


500 


reconstructed 


4.003 x 10 5 


3.704 x 10 5 


0.078 


1023 


822 


525 


error 


+2.6% 


+2.9% 


-0.25% 


+2.3% 


+2.8% 


+5.0% 



Table 15. The result of mass reconstruction from fitting the peak points for Example 2 
in Appendix A. The errors are calculated using rcconstr ucted-true ( except f or l g(Ai/A 2 ), 
which is reconstructed — true) and do not represent the statistical fluctuation. The statistical 
fluctuation are expected to be at the same order as the example in Section 2. 



Example 3 

This is also a simplified model with vector-like couplings. The mass spectrum and spin 
configuration are summarized in Table 16. The invariant mass distribution is shown in 
Fig. 24. The scatter plot and the peak-point plot are shown in Fig. 25. The result of 
mass reconstruction from fitting the peak points is shown in Table. 17. 





Y 


X 


N 


spin 


1/2 


1 


1/2 


mass[GeV] 


800 


300 


100 



Table 16. The summary of the mass spectrums and spin configurations of Example 3 in 
Appendix A. 
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Figure 24. The histogram of the invariant mass-squared of the two visible particles for 
Example 3 in Appendix A. In the massless limit (p± + P2) 2 = 2piP2- The vertical line 
indicates the value of AlA2 . 




Figure 25. The scatter plot (left) and the peak-point plot (right) of log {E\t/E2t) vs. 2pip 2 
for Example 3 in Appendix A. In both case the red curve is log (E\t / E2T) = log ( AlH ^^ lP2 ) 
with true Ai and A2 and the blue curve is the fitted one. 





AxfGcV 2 ] 


A 2 [GeV 2 ] 


log(Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


mjv[GeV] 


true 


8.000 x 10 4 


5.500 x 10 5 


-1.93 


800 


300 


100 


reconstructed 


8.480 x 10 4 


5.656 x 10 5 


-1.90 


815 


313 


115 


error 


+6.0% 


+2.8% 


+3.0% 


+1.8% 


+4.4% 


+15% 



Table 17. The result of mass reconstruction from fitting the peak points for Example 3 
in Appendix A. The errors are calculated using rcconstr ucted-true ( except f or l g(Ai/A 2 ), 
which is reconstructed — true) and do not represent the statistical fluctuation. The statistical 
fluctuation are expected to be at the same order as the example in Section 2. 
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Y 


X 


N 


particle 


2nd chargino 


left handed anti up squark 


1st neutralino 


mass[GeV 


1000 


770 


140.5 



Table 18. A summary of the decay chain in Appendix B.l. 



B Problematic Cases for Mass Determination 

In this Appendix we discuss cases where the mass determination method does not work 
even at the parton level. The two cases mentioned in sec. 2 are discussed in more details 
below. 

B.l The mother particle Y is polarized and decays anisotropically 

Our mass determination method relies on that the distribution of log *~!~o y cos ^ iy I is 

° 1+Py cos U2Y I Y 

symmetric around zero. This is obviously not the case if the decay of particle Y is not 
isotropic and preferentially emits particle 2 in the forward (or backward) direction. It 
can happen if the particle Y has nonzero spin and is polarized. Even if Y is polarized, 
sometimes a P or CP symmetry may still ensure the decay to be symmetric when 
all final states are included, as in the case of a Majorana neutralino, so some chiral 
structure of the process needs to be present. In the minimal supersymmetric standard 
model, the only particles that can have asymmetric decays are the charginos, so the 
study case that we consider is a decay chain initiated by a chargino. The particles in 
the decay chain are summarized in Table 18. The underlying model is SUSY LM2, with 
the mass of the heavy chargino adjusted by hand so that the decay chain can occur. 
To illustrate the polarization effect we restrict the heavy chargino to be produced only 
from the t-channel processes to obtain a larger polarization. 

The average value of the cosine of the transverse angle between particle 2 in the 
rest frame of Y and particle Y in the lab frame (cos^y \ Y ) is 0.124 (while (cos#iy | y ) = 
—0.033). The distribution of log i+^y cose^y \y a ^ er we DO °st eacn event to the frame 
where E 1 /E 2 = E 1t /E 2 t is shown in Fig. 26. The distribution has a bias and the 
average value is —0.136. The scatter plot and the peak-point plot (right) are shown 
in Fig. 27. The "peak points" are clearly shifted downwards. The slope of the fitted 
curve is smaller than the true slope, but the intercept is still close to the true value. 
The result of mass reconstruction shown in Table 19 exhibits large deviations. 

B.2 Soft visible particles and off-shell effects 

As mentioned in sec. 2, if the mass difference of two neighboring invisible particles in 
the decay chain is small, the off-shell effects and the "fake polarization" created by 
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Figure 26. The distribution of log | y (Left) and log i+^ToseZ W vs - invariant 

mass-squared (Right) of Appendix B.l. Each event is boosted longitudinally to the frame in 
which f = fe. The mean value of log ]ttZoZ \v is -°- 136 - 




Figure 27. The scatter plot (left) and the peak-point plot (right) of log (EytIE^t) vs. 2pip 2 
for Appendix B.l. In both case the red curve is log (E\t/E2t) = log ( Al +^ lP2 ) with true Ai 
and A2 and the blue curve is the fitted one. 





Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log(A!/A 2 ) 


m y [GeV] 


m x [GeV] 


m7v[GeV] 


true 


5.732 x 10 5 


4.071 x 10 5 


0.342 


1000 


770 


140.5 


reconstructed 


1.147 x 10 6 


8.092 x 10 5 


0.349 


1780 


1536 


1101 


error 


+100% 


+99% 


+0.67% 


+78% 


+99% 


+683% 



Table 19. The result of mass reconstruction from fitting the peak points for Ap- 
pendix B.l. The errors are calculated using rcc ° nstr uctcd-true ( excep t f or log(Ai/A 2 ), which 
is reconstructed — true) and do not represent the statistical fluctuation. The statistical fluc- 
tuation are expected to be at the same order as the example in Section 2. 



the Et cut may cause our mass determination method to fail. Here we examine these 
effects in more details. We first consider a decay chain with a small A 2 . The spectrum 
is summarized in Table 20. The on-shell approximation does not seem to work well for 
particle Y, which can be seen in Fig. 28. The distribution of 5 log ( Al + 2 ^ 2 ) (denned 
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Y 


X 


N 


particle 


heavy neutralino 


slepton 


lightest neutralino 


mass[GeV 


468 


441.9 


140.5 



B of events 



Table 20. A summary of the decay chain in Appendix B.2. 

tt of events 
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Figure 28. The histograms of yPy an< ^ y Px °^ Appendix B.2. The mass and width of Y 
are 468 GeV and 3.115 GeV. The mass and width of X are 441.9 GeV and 0.4076 GeV. 
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Figure 29. The scatter plot of <51og ( Al+ ^ lP2 ) vs. 2p 1 p 2 /^ z . <51og ( AlH g lP2 ) is defined as 
the effective value minus the true value, where the effective value for each event is calculated 
using py and p\ while the true value is calculated using m Y and m 2 x . The distribution is 
asymmetric, especially near 2p\p 2 / ^z' 2 = 1- This is due to a strong correlation between the 
maximum invariant mass and the effective values of Ai and A2. 



as the effective value minus the true value, where the effective value for each event 
is calculated using p\ and p x while the true value is calculated using my and m x ) 
vs. invariant mass-squared is shown in Fig. 29. One could see that it is asymmetric 
especially near 2pip 2 = A ^t 2 ■ This is due to a strong correlation between the maximum 
invariant mass and the effective values of A x and A 2 . 

In addition, particle 2 tends to be very soft in this example and the cut on Et 
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E T cut (GeV) 


(cos% y) (transverse) 


(cos0 2 y|y) (3-D) 





-0.025 


-0.039 


10 


0.041 


0.011 


20 


0.269 


0.176 


30 


0.439 


0.315 



Table 21. (cos#2v|y) (transverse and 3-D) for the spectrum in Appendix B.2 with different 
Et cuts. 




Figure 30. The scatter plot (left) and the peak-point plot (right) of log [EytI^t) vs. 2pij>2 
for Appendix B.2. In both case the red curve is log (E\t / E2t) = log ( Al +^^ 1P2 ) with true Ai 
and A2 and the blue curve is the fitted one. 



introduces a fake polarization for the particle Y decay. In this 10 GeV E T 

cut is applied to the visible particles and the resulting mean transverse angle between 
particle 2 in the rest frame of Y and particle Y in the lab frame (cos^y \ Y ) is 0.0406. 
To demonstrate this fake polarization induced by the Et cut, we generated a few 
more examples with the same spectrum but different Et cuts. The results are shown 
in Table 21. It is clear that when the visible particles are soft, a large Et cut will 
introduce a large fake polarization which could result in poor mass determination. 

The scatter plot and the peak-point plot are shown in Fig. 30. It is clear that the 
distribution is shifted downwards. Furthermore, because the slope is very small, the re- 
constructed values are very sensitive to fluctuations which brings additional difficulties 
to the reconstruction. The result of mass reconstruction is shown in Table 22. 

As a comparison, we look at the same case, but with the decay widths changed 
to very small values to force on-shell decays, and also with no Et cut. The result is 
shown in Fig. 31. It is clear that the bias presented in the previous case is removed. 
This shows that indeed the Et cut and the off-shell effects are the sources of the bias. 

The effects discussed above can also occur when Ai is very small. To explore this 
case we generated events according to the spectrum in Table 23. In addition, the width 
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Figure 31. The scatter plot (left) and the peak-point plot (right) of log (Eyr/Eyr) vs. 2pxP2 
of the case with very narrow width and no Ep cut in Appendix B.2. In both case the red 
curve is log [Ext/E^t) = log ( Al +^ lP2 ) with true Ai and A2 and the blue curve is the fitted 
one. The fit to the scatter plot works well. For the peak-point plot, the fluctuation is large 
due to a small slope. However, it is clear that the bias in Fig. 30 is not present in these two 
examples. 





Ax[GeV 2 ] 


A 2 [GeV 2 ] 


log (Ai/A 2 ) 


m y [GeV] 


m x [GeV] 


m N [GeV\ 


true 


1.755 x 10 5 


2.375 x 10 4 


2.00 


468 


441.9 


140.5 


reconstructed 


3.051 x 10 5 


4.478 x 10 4 


1.92 


827 


800 


579 


error 


+74% 


+89% 


-8.2% 


+77% 


+81% 


+312% 



Table 22. The result of mass reconstruction from fitting the peak points for Ap- 
pendix B.2. The errors are calculated using rcconstmcted-truo ( except for i og (Ai/A 2 ), which 
is reconstructed — true) and do not represent the statistical fluctuation. The statistical fluc- 
tuation is large due to a small slope. 



of particle X is changed to about 4 GeV to enhance the off-shell effects. The scatter 
plot and the peak-point plot are shown in Fig. 32. The result of mass reconstruction 
is shown in Table 24. Although the value of A x is very small, the result seems to 
be acceptable. This is due to the following reasons: 1) although in the rest frame of 
particle X the available energy for particle 1 is small, in the lab frame particle X is 
more boosted on average than particle Y so that particle 1 is less likely to be very soft; 
2) the Et of particle 1 is less correlated with the relative direction of particle Y because 
it does not directly come from the Y decay; 3) the slope is large and the fluctuation is 
small compared with the previous case. 
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